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ABSTRACT 



The multi-loop Lead Anfle servo loop is utilized as a analog computing 
loop that mechanizes the range-bearing torpedo fire control equation. 

This thesis applies several new servo analvsis techniques, developes 
system component transfer functions, analyzes individual loops, and 
examines qualitative stability and response of the overall system. 

No response or error criteria was specified. Based on the author' 3 
experience and advice from submariners a minimum error criteria is 
4 derived. This thesis preposes several modifications of the origional 
system in order to improve performance, points out several difficiencies, 
and recommends areas of further study. 

The writer wishes to express his appreciation for the courtesies 
and cooperation extended to him by the personnel he had the pleasure 
of working with at the Librascope Division of General Precision, 
Incorporated at Glendale, California. 

The author also expresses his appreciation for the advice and encourage- 
ment received from Lt. Charles Poliak and especially Dr. Marvin Pastel, 
faculty advisor. 
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Definition 



Relative Target Bearing 



Missile Azimuth 



Torpedo Lead Angle 



Error in Bb6 
Relative Tube Bearing 



Torpedo Angle 



When used alone, B denotes 
the relative bearing to the 
present position of the target 
and is measured in the Master 
ovn ship horizon system. 

The difference between 
Relative Target Bearing and 
Relative Tube Bearing, mea- 
sured in the horizontal plane 
of the own ship horizon system. 

The angle formed by the 
intersection of the line from 
own ship’s present position 
to target's present position 
with the line from torpedo's 
pseudo or real position to 
target's future position. 
Positive angles are measured 
clockwise from the line connect- 
ing own ship and target, i.e., 
the present relative bearing 
line. Prior to event ii, and 
for all Preset problems, 
point li is the pseudo position. 
For Wire Guided problems, 
after event li, the dead reckoned 
torpedo position is used. 



Relative bearing of the torpedo 
tube measured in an own ship 
horizon system. 

Angle between the horizontal 
projection of the missile or 
torpedo speed vector and the 
present bearing line from own 
ship to the missile or torpedo, 
measured in the horizontal 
plane of a missile or torpedo 
horizon system. 
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Target Angle 


Angle between the target speed 
vector and the present bearing 
line from own ship to target, 
measured in the target horizon 
system. 


Co 


Own Ship Course 


Course of own ship measured 
in. the Master system. 


Cm 


Torpedo Course 


Course of the torpedo or 
missile measured in the 
horizontal plane of a north 
oriented horizon system at 
The torpedo or missile. 


j (Cm) 


Computed Correction 
to Torpedo Course 




Cmb 


Torpedo Co .rse at Fvent 

u 


The course of the torpedo 
measured in the horizontal 
plane of a north oriented 
horizon system at point b. 


Cmob 


Relative Torpedo Course 
at Event b 


The course of the torpedo 
measured in the horizontal 
plane of an own ship oriented 
horizon system at point b. 


Ct 


Target Course 


Course of the target measured 
in the target north oriented 
horizon system. 


D3y 


Bearing Rate 


The time rate of change of 
true target bearing. 


DCo 


Own Ship 'Turning Rate 


The rate of change of own 
ship course with respect 
to real time. 


DCt 


Target Turning Rate 


The rate of change of target 
course with respect to time. 


DKho 


Own Ship Speed 


The time rate of change of 



the distance traversed by 
own ship, measured in the 
horizontal plane of the 



Master own ship horizon 
system. 
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Symbol 

DFht Target Speed The time rate of change of 

the distance traversed by 
the target, measured in the 
horizontal plane of the target 
horizon system. 

f Frequency, cycles per 

second 

Fc Closed loop transfer 

function 

Fo Open loop transfer 

function 



G 


Gyro Angle 


Angle between vertical 
plane thru line of fire and 
vertical plane thru desired 
missile centerline measured 
in horizontal plane. (When 
used alone only.) 


e(G) 


Error in Gyro Angle 




sq(G) 


Angular Spread 


A manual correction to gyro 
angle to initiate an angular 
spread of torpedo paths. 


G6 


Gyro Angle at Event 6 


The torpedo gyro angle at 
event 6. 


Hm 


Lenpth of Torpedo Path 


Actual distance traveled by 
torpedo from own ship to 
target. 


c(Hm) 


Computed Length of 
Torpedo Path 




ec(Hm) 


Error in Computed Length 
of Torpedo Path 




q(Htn) 


Run to Burst Command 


Manual correction to torpedo 
path length. 


sj (Hm) 


Correction to Torpedo 
Path Length 


A bias added to torpedo path 
length to compensate for 
torpedo acceleration to 
required speed. 
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Name 


Definition 


sq(Hm) 


Enabling Run Offset 


Correction used to shift 
enabling point along the 
torpedo track. 


Hm26 


Total Run to Burst 


The total theoretical length 
of the torpedo path from event 
2 to event 6. 

(Obsolete symbol Hmt) 


Hm6 


Di stance from Torpedo to 
Point 6 


The straight line distance 
from the torpedo's present 
position to point 6. 


Hm7 


Remaining Run to Burst 


The total distance the torpedo 
must run from its present^,,,.- 
position to its bursting: position 
(point 7). 


Hvm 


Torpedo Run Depth 


Vertical distance the torpedo 
is set to run below the surface 
of the water. 


Hvmg 


Depth Difference 


Vertical distance between the 
torpedo tube and the missile 
or torpedo. (Hvm - Hvg) 


K 27 


Distance from event 
2 to event 7 




H67 

A 


Linear Spread Along Tar- 
get Track 


The distance along the target track 
between the future target position 
(point 6) and the advance point (point7) 

r 

✓ 


3 

Ly2 


Latitude of Firing 


The latitude of own ship at the time 
of fire. 


Lyo 


Latitude of Own Ship 




Lyp 


Latitude of Proofing 


The latitude at which the torpedo 
was proofed. 


Lyt 


Latitude of Target 




k 


One Thousand ohms 




K 


Gain 





X 



Name 



Definition 



Symbol 



Mho 


Own Ship Horizontal 




Motion 



Pdn 


Athwart ship Parallax 




Displacement 



Pdo 


Centerline Parallax 
Displacement 


r ^ , r 2 . * . . 


R~ots of equation 


R 


Present Slant Range 


j(R) 


Computed Correction to 
Range 


Rh 


Horizontal Range 



c(Rh) 


Computed Horizontal Range 


j(Rh) 


Horizontal Range Correction 


Rh2 


Horizontal Range to Target 




at Event 2 



The horizontal distance made good 
by own ship during a given inter- 
val of tine. 

The component, in the deck re- 
ference plane and perpendicular 
to own ship centerline, of the total 
distance from the origin of the 
Master system to the origin of the 
weapon system. (The center of the 
torpedo gyro or the missile platform.' 
Pdn is measured in the Master own 
ship deck system. 

The component, along own ship center- 
line, of the total distance from the 
origin of the Master system to the 
origin of the weapon system. (The 
center of the torpedo gyro, or the 
missile platform.) Pdo is measured 
in the Master own ship deck system. 



R, when used alone, is the 
present distance between 
own ship and target, mea- 
sured along a straight line 
in the Master system. 



The present horizontal range 
from own ship to the target's 
present position as seen by 
radar or as computed from 
sonar data, measured in the 
Master system. 



Computed correction to 
horizontal range. 

The range, measured in 
the Master system, from 
own ship to target at time 
of fire (event 2). 
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Definition 



Py mbol 

Rh3 



t 

c< 

fi 

G 

8 

Q 

Q+ 

Go 

Tf 

UJ 

7 

Mf 

T26 



T 

l x 

Urn 



j(U») 



Name 

Horizontal Rr-nge to Reach 
Point 

Time 

TVs 

Base of Natural Logarithm 
Frror 

Angular displacement, '"radians 

System Input 

System Output 

Time of Flight 

Angular Velocity, radians 
per second 

Time Constant 

Microfarads 

Time of Travel 

Time to Event x in 
Fig. 1. 

Torpedo or Missile Speed 



Correction to Torpedo 
Speed 



Tie range, measured in 
the Master system, from 
ovm ship present rosition 
to the torpedo reach point, 
point 3. 



The total time from time 
of fire (event2) to time of 
torpedo burst (event 6). 



The speed of the torpedo or 
missile with respect to the 
medium through which it is 
moving. The time rate of 
change of the distance over 
which it is moving. 

A computed correction to 
torpedo speed. 
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Symbol 



Name 



Definition 



Ym 



Torpedo Turn Radius 



Target Half-Length 



The mean radius of curvature 
of the path taker, by the torpedo 
(essentially circular) in 
traveling from the reach point 
(point 3) to the gyro involute 
point (point U) . 



Distance from the stern of the 
target to the middle of the target 

(MOT). 
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I . SUMMARY 



An analysis of individual loops was conducted and a representative 
loop was selected and presented herein. The individual loop as 
designed will operate satisfactorily; however, using optimum settings 
of input and feedback resistances determined in section 5-1 the response 
can be improved 37.5 per cent over the best response obtained by trial 
and error, using settling time as a criteria. A further imorovement 
was obtained by eliminating the integral-damping network; however, 
considerable noise was incountered which would not show up in the output 
but would cause gear wear in the chain near the servo motor. By chang- 
ing capacitors C-3 and C-U in figure lit to .25 microfarads instead 
of the presently installed 1 microfarad, the same resnonse could be 
obtained without the undesirable oerformance at steady state. The over- 
all imorovement in response using the resistances determined and the 
change in the caoacitance was 5b. 2 per cent. 

The overall system wa3 considered but because of the non-linear nature 
of the resolver in the circuit exact response cculd not be predicted in 
the time allowed for analysis. Linearizing assumptions were made and a 
stability analysis was performed. It was found that the system was high- 
ly stable with a gain margin of 2b db and a phase margin of 72 degrees. 
During the analysis it was discovered that the system was type zero 
which has an infinite steady state error to a ramp input. Within the 
area of expected performance of the system this error did not appear to 
be significant; however, if large ramp inputs are incountered or long 
C'mputing times the error can become significant. The system is easily 
compensated to become type one by an additional amplifier-motor loop. 



1 



If it is desired to use the system as designed it appears that by 
careful scaling that the error can be reduced significantly without 
further compensation. 
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II. NTR'' DUCT I ON 



2. 1 The Gyro Lead Anple Conputinp Loop 

The pyro lead anple corputinp circuit functions as a part of the 
Attack Director, Mk ?5* It is corsiiered advantapeous to familiar ~e 
the reader with the overall system prior to discussinp the specific 
circuit analyzed. 

The Attack Director, 75 performs as a torpedo data computer and s : t- 
icn keeper. Cortinr :s inruts of ran re and bearinp information a r - rr- 

from various selected data rovers. 

* 

As a oosition keeper the Attack Director renerafes ranpe and Dear ns of 
tarpet positi n for transmission to sonar v»a the Analyzer Console, * 
fro^ the most recent values of tarpet course, soeed, rarpe, and bea 1 " r 
r peneratos ranpe and r 'earinp of apparent tarret position for direct 
trarsmission to sonar, from the nest recent val jes of automatically r 
manually inserted inputs of tarpet course and speed. 

As an anp e solver the Attach Director computes pvro arrle, orat , ‘ s , 
ar d t do r nninp erth f^r electrically set torredoe^.. 

p quired information r or the -* rmrol of wire ru led torpenceo 

9 

-regions to torpedo c' ur se and rm to turst, and a ^ertir; + 
f r >. ’rv. i s acc r i si f.*d r o : * c + r. 

" ! of ovr shir •* n ‘ r shif a^e rientateil r ; ' ♦ . j- 

np f he v easin ' lire to a t ^ et # ^ rredo inputs arr r .r' 
o' n ^ior to ti^e of' fir* 'V r m ruined and pres t t ,r r •••! * /. 
mr re rcsition is nspla- r ed a*’ter tt^e of fire for th* vire ’ r 
~ub : ns. 



n 



The Attack Director, Mk 75 functions in both automatic and display modes. 



The specific loop which is the subject of this analvsis is a multi-loop, 
electro-mechanical servo-mechanism. The loop is composed of eipht individ- 
ual loops, which are interdependent in the operation of the Attack Director, 
Mk 75. Listed below are the element loops of the compensation of lead 
anfle loop and their phvsical location in the director: 

1. Gyro Anfle - located in the Angle Solver Section 

2 . Enabling Run - located in the Angle Solver Section 

3. Torpedo Lead Angle - located in the Position Keeper Display 

Section * 

Ii. Relative Target Bearing - located in the Position Keeper 
Display Section of the Attack Director 

5. Relative Tube Bearing - located in the Ballistic Computer 
Section 

6. Weapon Time to Purst - located in the Ballistic Computer 
Section of the Attack Director 

7 . Target Speed - located in the Position Keeper Display 
Section 

8 . Target Angle - located in the Position Keeper Displav 



The commutation of the lead angle loop is a sequential computation loop. 
This loop has four possible mechanical inputs and five possible electrical 
inputs. The outputs from this loop are three dial readings. These 
readings are valid onlv when the loop has reached a steady state condit- 
ion. 

In order to explain the operation of the computation of lead angle loop 
several assumptions were made. Initial conditions are such that Gyro 



Section 



Angle [C] , Torpedo Lead Angle [8t> Lj , Enabling 




, and 



b 



Weapon Time to Burst Kib are initially at some particular value. Another 
condition is that electric cower to all loops is not applied until such 
time as the computation begins. 

Another assumption made is that all electrical inputs are constants with 
the exception of the following variables. 

1. Ym' Sin Bb6 

2. Cos Bts (Dmht -K1 . T26) 

1. Ym' Cos Bb6 

L. Sin Bts (Dmht . K1 . T26) 

With the above in mind let us consider the solution to the problem: 

The firing of a preset torpedo reouires solving for a predetermined point 
of impact between torpedo and the target, noted as point 6. The torpedo 
innut requirements are: run to the enabling point, c(Hm); the gyro angle 

order, G; run depth, Hvm; and a 26 Volt mark signal when 0°<T<180^ for 
automatic horizontal search. 

Computations are made to resolve into mutually perpendicular components 
the final vector (leg) of the torpedo run to the target. This is the 
straight line connecting event h with event 6. The components are JETY, 
which is parallel to the line of sight (LOS) and 2X, which is perpendicular 
to it, (reference figure 1). 

The ouantities Z X and £Y are predicted by two isolation amplifiers summing 
the applicable inputs. These two amplifiers are the start of a seouential 
computation. The amplifier inputs involve the parallax quantities of the 
tube, Pdo and Pdn, the reach of the torpedo Rh3, the torpedo turn radius 
Ym' , the relative Target Bearing B, corrected for Relative Tube Bearing 
Bp, the Torpedo Lead Angle Bb6, the Target Angle Bts, target speed times 
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the total run time of the torpedo (DMht . T?6), and a linear spread H67, 
quantity if cresent. These quantities are related to £X and £Y by the 
following equations. 



X = sin B (Pdo+ Rh3 cos Bp - Ym' sin Bp) 

- cos B (Pdn + Rh3 sin Bp - Ym' cos Bg) 

+Ym' cos Bb6 + sin Bts (DMht . . T26 + Ho7) 

Y = Rh - cos 3(Pdo + Rh3 cos Bp - Ym' cos Bg) 

- sin B (Pdn + Rh3 sin Bp - Ym' cos Bg) 

- Ym' sin Bb6 - cos Bts (DMht . K x . T 26 + H6?) 



ver. A servo loop positions the mechanical innut of the resolver until 
one of the rotor outputs approaches null. The output of the other rotor 
is the vector sum of £X and £Y Hm6. The shaft position of the servo loop 



A perturbation, which includes a compensation for drift and gyro uncag- 
ing time, is summed with the null signal e(Bb6) from the resolver in an 
isolation a"plifier. The generation of this perturbation will be explained 
later. The output of this isolation amplifier is sent to a variable gain 
amplifier, VGA, which divides this signal by the other output of the lead 
angle resolver. 

The output of this amplifier is the error in Bb6. 



The cuanti ties £X and 2-Y are inputs to the two stator windings of a resol- 



is an angle equal to the tan " ' ~~ « Bb6, the torpedo lead ar.ple (refer- 

y 

ence figure 1). 
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This error is computed by the relation 9 “pwhere 0 is th^ error in Bb6, 
s is an error in the difference in length of a particular component of 
£ X and £ Y and the perturbation mentioned above, and r is the vector 
sum of 2. X and X Y, Hm6. 

The vector sum of IS X and £T is the length of torpedo travel from point 
I* to point 6. This signal is sent another isolation amplifier, whose 
function will be explained later. The shaft position of torpedo lead 
angle (Bb6) is displayed on the outer dial of the "Target Group" and "Own 
Ship Group" dial display in the Position Keeper Display Section. 

The Bb6 shaft position is added to '.he target angle (Bts) by a mechanical 
differential (Bb6 + Bts * I). The output of this differential positions 
a limit cam, which operates a micro switch. The normally closed Dosition 
of the microswitch gives a right turn mark signal to the torpedo for 
automatic horizontal search, NC, 0<I<180 . The normally open position 
provides no mark signal and therefore, aleft turn is in order. This turn 
signal is available as a oreset torpedo function that is activated at 
the enabling point. 

The Bb6 loop also positions a two speed control transmitter synchro sy- 
stem. These CX's send an electrical analog signal proportional to the 
Torpedo Lead Angle, Bb6, to a two apeed control differential transmitter 
svnchro system that is raechaically referenced with the Relative Target 
Bearing, B. These CDX's add the two functions and yield Relative Torpedo 
Course, Cmoln The CDX's then send this signal to the Own Ship's course 
module. The Cmoli signal al3o references a two apeed control differential 
transmitter synchro system that is mechanically referenced with Relative 
Tube Bearing, Bg. The output of the relative tube bearing CDX's reference 
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a two speed CDX svstem on the Angular Spread module sq (G). This ~odule 
is located on the firing panel. The output of these CDX's reference a 
two speed control transformer svstem on the Gyro Angle Order module, G. 

These CT's are the feedback element for the loop. The equation involved 

here is: [G = G04- + ^(CrJj . 

Each of the terms on the right side of the above equation for G have been 
added by the use of two speed svnchro systems (CDX). This signal refer- 
ences the gyro angle module feedback CT's. These CT's have a shaft position 
rropcrtional to Gyro Angle, G, and an output, which through a servo amp- 
lifier and motor, positions the gvpo angle order shaft. 

The pyro angle order gear train positions the wiper on a DC potentiometer 
(rot). This G pot is referenced with TTtimes torpedo right turn radius, 

77- y ^ and 77* times torpedo left turn radius, 77* Y m ^ . The center 

tap of the pot is grounded. The voltage from the not is positive for both 
ripht and left turns, but it is not of the same magnitude. The voltage 
output at the wiper position of this pot is then proportional to the wiper 
position, a function of G, and the reference voltage. The analog signal 

is then the product of these two inruts; hence the arc distance traveled 

77 " 

about the gvro angle, Ym' . G . This signal is sent to an isolation 
amplifier, figure 2. Other inputs to this amplifier are: The distance 

to the reach point, Rh 3 , from the ballistic constants network; the vector 
sum, X sin Bb6 + Y cos Bb6 ■ Hm6j , from the lead angle resolver 
previously discussed; and a correction for torpedo run depth Hvm, and 
own ship depth Hvo, Hvmg^ , from the ballistic constants network. 

The output of this amplifier is proportional to Hm26. This voltage is 
sent to an impedance matching transformer in the input network of a servo 
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ISOLATION AMPLIFIER SCHEMATIC 



amplifier. The output of this network is the output of the feedback 
vernistat Hm26 on the Enabling Run servo module. The amplifier drives 
the servo motor until the two inputs are equal in magnitude and 180 out 
of phase. The motor has then driven out a mechanical analog position 
proportional to Hm26. This shaft position may be added to or subtracted 
from manually by a hand crank input labeled Enabling Run Offset, sq(Hm). 
This is accomplished by the use of a mechanical differential. The amount 
of offset is recorded on a counter for display. 

The Enabling Run module visually displays computed Enabling Run, c(Hm), 
which is the summation of Hm26 and, the offset sq(Hm). One vernistat on 
this loop sends a voltage proportional to the computed Enabling Run c(Hm), 
to the Run to Burst module as a reference source for a vernistat feedback 
loop. This particular total run to burst, [Hm26 w/g - NR ■ c(Hra) w/g + 
q(Hm)J , is used only in the wire guided solution and will not be discussed 
here. 

A separate cam for each weapon on the computed Enabling Run shaft, c(Hm) 
figure 3, sends a within limit signal to the gyro angle limit cam. If 
the gyro limit cams are within limits, then they send a signal to the 
depth stratum limit cams. The run depth cam, then sends a sipnal to the 
within limit circuitry in the weapon monitor unit. 

The Enabling Run module position six one speed control transformers that 
send an analog error signal, ec(Hm), to a servo motor in the torpedo. 

There is a control transmitter, mechanically linked to a servo motor in 
each torpedo. This CX feeds a follow-up signal oc(Hm) back to the CT's 
in the Enabling Run module to comolete the servo loop. 

The second vernistat in the Enabling Run module is referenced with the 



11 



/Vo r£ 3 



J. TH£X£ AX£ //V O/v/O c Al C T'S 

fO R. £ A c >V 7 ~ L/ 3 

2 . A r/zexe AAf z+o/y/o vaz. 

fax £AC/Z TVP£ of >*/£APQM. 





TO ^ 6£ \* / rP//v L />m / t~ 

l/qHt or* F/Ar cozv. S£c. 
7 H*U BAH. 3 £L. SW- 



4 2 8 v DC 



(S/rCVi.- SuvLYf>)rZL-<.(Hm)/* J 



y/GOA? 3 

yyyy&t/A' g /&u/r 



it 



+ D CO . To • c ( Hf* ) 




difference between the sine of firing latitude and the sine of proofing 
latitude times half the product of the enabling run time and target turning 
rate, this minus the signal for gyro uncaging time, times own ship turn- 
ing rate. These two quantities are multiolied by enabling run by the use 
of the vernistat. 



r / 

T2L { S//y l ip ~ L y 2 ,~ oc °' 

The output of this vernistat is sent as the perturbation to the lead am le 
servo loop, figure b, previously discussed. 



The analog voltage, c(Hm), from the Enabling Run Vernistat, is sent to 
the total run to burst module, figure 5, and to an isolation amplifier, 
which sums computed Enabling Run, c(Hm), and a torpedo path length correct- 
ion, Hm corr 



Hm 

corr . 



|sg(Hm) ± Hvmg (Hm) + Hvmf 




(Note: The Hvmp^p^) t ' erm USec * only with the Mk 37 Mod 0, high speed 

torpedoes.) The output of this amplifier is total torpedo run, (-) Hm27. 
This signal is sent to an impedance matching transformer in the in'-ut 
network of the isolation amplifier of the Weapon Time to Burst servo systen. 
A vernistat on the Weapon Time to Burst module is the other inrut to the 
amplifier. This vernistat is referenced with torpedo or missile speed. Urn 1 , 
from an isolation amplifier having 3 ballistic constants f(Hmv), Ur, i(Um) 
for inputs. The f(Hvm) signal, a function of depth, is a correction to the 
torpedo or missile speed, Um. This compensation is generated on the Run 
Depth module by a pot positioned proportionally to Run Depth, Hvm, and 
referenced from the ballistic constants network with a function f(Hvm). 

The voltage output of the T26 conputing vernistat is proportional to 
torpedo speed Um' times weapon time to burst T26. When the product of 
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of speed and time (ecual to total run distar ce Hmp), is compared with 

( - )H27 and passed throuph an isolation amplifier, it is an error signal, 

e(T26), used to position the Weapon Time to Burst, T26, servo system 

which displavs T26. A vernistat referenced with the difference between 

the sine of proofing latitude and the sine of firing latitude times half 

the torpedo turning rate and positioned by the T26 servo system has an 

analog voltage output equal to the rroduct of weapon time to burst and the 

fDOm . , .] 

compensation for firing latitude . T26 (sin Lyp - sin Lv2)J . This 

signal is then summed at the input of an isolation amplifier with a con- 
tribution from the ballistic constant network. ( -DCo . Td), where DCo 
is the Own Ship Turninr Rate and *Td is the Gvro Uncaging Time. The out- 
out of the isolation amplifier is used to reference the vernistat on the 
Enabling Run nodule previously discussed. 

Another vernistat in the Weapon Time to Burst module sends a signal 
proportional to T26 to reference a vernistat on the Target Speed module. 

This vernistat' s output is a multiple of Target Speed, DMhtj Weapon Time 
to Burst, T26 and K]_, to yield the distance the target travels from time 
of fire to time of impact. This signal is sent to an isolation amplifier, 
figure 2, that sums this distance and linear spread quantity H67, if 

present. The linear spread inDut is located in the firing nanel. Linear 

spread is never used concurrently with Angular Spread auantities, sq(G), 
previously discussed. The output of this amplifier, the total distance 
travelled, is [(Tf . DMht . %) + H 67 ] or [(T26 . DMht . K x ) + H67 J , 

and references a resolver in the Target Angle servo system. The output 

of this resolver is the sine and cosine of Target Angle, Bts, times the 
reference signal. These two signals become part of the inputs to the 
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two isolation amplifiers which vield F X and £ Y. 

Generation of one of the other inputs to each of these amplifiers bepins 
with the output of two isolation amplifiers. The input to these amplifi- 
ers is dependent uoon the value of the Gyro Angle, G. Three conditions 
can be considered: 

(1) For 170°^ G ;>10° and l^O^ G >190° the inputs into the turn- 
radius-isolation amplifier is just Ym' and the inout into the 
reach-point-isolation amplifier is Rh3 + ^^^Rhl w here Hvmgj^^ 
is the correction to Rh3 for depth difference, (Hvo - Ph - Hvm). 
Under these conditions and £Y are computed as: 

Z X » [pdo sin B - Pdn cos B + (Rh3 + Hvmgj^) sin (B - Bp) -ym' 
cos (B - Bp) + Ym' cos Bb6 + (DMht . T26 . K x + H67) sin Bts J 

(2) For 350°< G < 10° the incut to the pvro-turr-radius-isolation 

amplifier is zero; however, the reach-noint-isolation amplifier 
's now Rh3 + + (Ym 1 . G . — ^/qq ~ ^ af ^ition 

to this change, additional circuitry removes (-Ym" sin Bbfc) 
and (Ym' cos Bb6) from the £X and £Y isolation amplifiers. 

Hence £ X and £ Y are commuted as: 

£ X = Pdo sin B - Pdn cos B + (Rh3 + Ym' . G . 1 ) 

s&O* 

£ Y = Rh - Pdo cos B + Pdn sin B + (Rh3 + Mvmp PVl -j + Ym' G. — ) 

ru > 0 

cos (B - Bp) + (D>'ht . T26 . + H67) cos Bts 

(3) For Y'C a <. G < 190° the innut to the pyro- turn-radius-isolation 
amplifier is 2 Ym' or twice the value used in case 1 above. The 
reach-noint-isolation amplifier is the same as in case 1 and as 
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oefore, additional circuitry removes (Ym' sin Bb6) and Ym' cos 
Bb6) from the SIX and 51 Y isolation amplifiers. 

Furthermore, it should be noted that the quantity Ym' used in 
(1) is actually Ym^ or Ym r , depending on whether the condition 
360^C^180 exists or not. 

In case (2) the same thing is true, however, its determination 
is due to the position of a not positioned in the Gyro Angle 

module, figure 6. 

In case (3) Ym' implies 2Ym'^ or 2Ym' r , depending on how the Gyro 
Angle develops. 

These signals described above for case one, two, and three, reference 
the stators on a resolver on the Relative Tube Bearing nodule. This 
module has two positions. One position is for the bearing of the 
port tube and other position is for the bearing of the starboard tube. 
The incut signals come from the own ship section where the.e two 
bearings are stabilized. The distance to the Peach Point (point 3) 
and gyro-turn radius are measured along and associated with the tube 
bearing module references these ouantities to own ship's centerline, 
(Own Ship's Course, Co.). 

To these quantities are added the appropriate parallax components, 

Pdo, Pdn, of the ship by summing the applicable input in two isolation 
amplifiers. The output of these two amplifiers then reference a re- 
solver on the Relative Target Bearing module. This resolver rotates 
the input quantities to distances along and across the bearing 
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line (LCo). These quantities are then sent to the summing amolif ; ers 
mentioned before and v ; ' eld Z X and Z Y. 



2.2 Objectives of Analysis 

The analvsis was undertaken to determine the response characteristics 
of the individual loops, the overall loop, check the system for stab- 
ility, recommend potentiometer settings for optimum response and if 
found necessary recommend compensation networks. 

2.3 Period of Time Allowed for Analysis 

Work on this thesis was accomolished during eipht weeks of the summer 
Deriod between the second and third year of the Ordance Engineering 
Guided Missiles Curriculum ending in December of i 960 . 
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III. Review or Introduction to Some Recent Developements in Servo 
Analysis Techniques 

3.1 Determinantal Method of Analysi s^ 

The theory of simple single loop control systems is well developed, and 
the block diagrams of such systems are not complicated. On the other 
hand, the block diagrams and theory of multi-loop, multi-input, multi- 
output, and multi-coupled control systems is not very well developed 
because of the complexity of the problem involved. The block diagrams 
of such systems are complicated and difficult to intepret. One approach 
to the solution of such a problem is to formulate and solve a set of 
simultanious equations for the system. This can best be accomplished 
by arranging such systems into a standard block diagram. Arranging the 
system in this manner greatly facilitates the writing of the equations. 

Any system which has summing points in the main transmission path can 
be arranged in a block diagram by suitable block diagram manipulation 
as shown in figure 7. Writing the equations at the various nodes: 

At node a 

a = R a * (G ab* G Aa) a “ (°bc # G Ba^ b " (G cd* G Ca) c 
A more convenient notation is: 

a • » a - <°aa> a ' <°ba> b * < G ca><= 

At node b 

» ■ *b - «W b ♦ (°ab> a ‘ < G cd )c 

^CHU, Y., A Generalized Theory of Linear Multi-Loop Automatic 
Control Systems , Doctoral Dissertation, MIT, 19^3 



% 
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7 5 7 -^/^ OCA. 





at node d 



d = (~ ) d + fa +(G a j)4 + d jJ c 



Rearranging to bring out the symmetry 



(j 1 G&a) CL +\.Gt)<x) b ~t \.G ca.) C +{£<ia)d - 

(- GclU a + (i + Gkb)b + (GcbJc + LCdbJtl Pb 

(- Gac) CL - (Gkc)b +(i+Gcc)C + (*jJc ~tc 

C- G«dJ 4 - ^ojjd ’ (G^i)c i-(l^0dfijd ^ td 



solving by Cramer’s 


Rule 
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For simplicity let us assume the inputs 

- szUj ; x * ^ ^ 

Then the output will be 0 ~ (.(?<* c ^ d 
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It is noted that A is equal to the characteristic equation of the sys- 
tern. For an example, take the four node system in figure 8. 
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equation of the system. 
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Now putting the transfer functions in the form of numerator divided by 
denominator. 



£ - 



k, /A z /As //y 

is i. /s * V v 



X, //- /A //. *v- 
L * ^ t i ' z ’ 



//z ^ ~ 



X?.: ^ * ✓ Xj - 



Which when simplified becomes: 






X? /Kt /fry / fry 

//'v XV XV _ 74 X/ A/* A/ I A/y xy 



/V* A/j A/* ASy 



jm 

The roots of the transfer function are the zeros of the characteristic 
equation. 



Oi A» £s- X *' A/j. Afj /Xy i)jr A /Vt, sjj /PV As r ~ O 



Now if we desire to determine the effect of the feedback path G5, we can 
rearrange the characteristic equation in the fora 

Oy /XU £>± XV a /c, A/x A/j A/ V j A A/z. A/j A/ v AAs ~ ° 

The quantity in the brackets can be equated to zero and by root locus 
techniques the roots can be determined. The characteristic equation will 
then be fis [ {6 +A,)(s + ^LxJ v- 5 * AsJfat ^ ^ ^ 
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The final root locus equation will then be: 



JY i //> A '■/ A/y 



- 4s 5" (S 



Ps(j+n.)lS-> n-^jiS +a,j srA y l t +■ r c w 

The transfer function will then be: 



- A 






£ x il£l Idlll l±L 
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3.2 Determination of the Composite Open Loop Transfer Function and 
Velocity Constant for ;4ulti-Loop Systems 
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In terms of numerator and denominator of the closed loop transfer function 

r. — — — 

£<. - /V» 

Now applying the classical definition of the velocity constant 



A* = 



Am. S A* 

S - 9 - o 

AA 5 /K. 



<^ c - 



It is convenient to determine the open loop transfer function in this 
manner and then apply root locus techniques. 
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After the open loop function has been determined, it is quiet easy to 
determine the type and order of the system. It will have the form: 

s _ A (S> (S * ijl • 

S «[ ( • ■ ■] 

The denominator of the open loop function is unfactoredj the S° can be 
factored out prior to performing the root locus or K versus real roots 
thus the type and number of the system can be determined. If the desired 
type is not obtained, this is a criteria for selection of suitable comp- 
ensation. In any event it avoids the necessity of finding the roots of 
the system to determine the type number. 



It is often necessary to perform a root locus or K versus real root plot 
to determine the poles of the open loop transfer function. As an example, 



consider the example of section 3 « 1 » The closed loop transfer function 



is 



/c = 
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The denominator then can be factored. 

Vi Os Vy Vs v /Vi /Vj /V* /Vj - v? 
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The roots are then determined and the open loop transfer then becomes 



/T - A 4 t J r ■A 'i ^3 A/y 

Lf/ AiJ 

The closed loop roots are then determined in the usual manner. 

3.3 Root Locus Gain versus Real Root s? 

The usual method of determining the roots of a closed loop function is to 
apply the familiar root locus principals. A somewhat similiar technique 
which embodies some of the same principals is the method of plotting 
K (root locus gain) versus the real roots. 

For illustration, let us consider this unity feedback loop shown below. 



- jF.- % j L 

The closed loop transfer function is 

r ^ OX =• 

/^/; 

The roots of Fc are where 0* A/s — o 0 r 

The open loop transfer function has this form* 

A/a_ _ * (sy2, j vS* £ *.) • 

0* S'* vi trtj 

2. Wheeler, R. C. H., Unpublished Class Notes 
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Now plotting the poles and zeros on the S - plane: 




By the Evan's method the closed loop roots are determined by summing angles 
to 180° and plotting the locus of the roots. Then knowing the gain associated 
with the system, a cut and try process is performed to locate the roots, by the 

relationship K~ —£± . This process is necessary along the real axis where 

Hi: 

root locus exists. If these values along the axis are plotted the following 
curves w ill result . All the curves shown do not result from every open loop 
transfer function but are shown as an example of all possible configeration. 
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Fig. 10 PLOT OF ROOT LOCUS GAIN 
VERSUS REAL ROOTS 

It can be shown that the maximum between and f * > is the "breakaway point" 
and the minimum between Z, and Z A is the entrance point in the Evan's Root 
Locus Method. 

If the system gain is given or can be chosen, a horizontal line drawn which 
intersects these curves yields the real roots of the closed loop function. 
If it does not intersect as between and // , complex roots exist. Also 
between Z 3 and /i , one real root and two complex roots exist with the 
gain chosen as shown. 
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This method works equally well with complex poles or zeros in the open loop 
function by using the same relationship for gain k — ~^L±_ , but one 

must remember that the distances to these complex quantities to the point 
at which the gain is being computed must be included. 

From the foregoing discussion, it is seen that all the real roots of the 
closed loop can be determined without plotting the root locus. 

It is a fact that the denominator of the closed loop function which is the 
characteristic equation of the system has the following form and possesses 
the below listed qualities: 

AX* t bx n ~' D = o 

A - 1 

B ■ Sum of the roots of the closed loop function 
C - Product of the roots of the closed loop function 

Knowing the real roots from the plot of "K versus Real Roots" and using the 
relationship A. + A*. + l j (oC +j c ) + ^ - J uj^J = £ 

+Aj + 3.4. * y 

we can now solve for the real part of the complex pair. Knowing this we can 
use the second relationship 

/X, n * x , i-*- tj 3 0 

which becomes 

X, /u a* ( «. * y- - U 

Since Z x /£ i AM 0 are known, we can solve for Wc thus all roots are 

known. 
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This method provides a criteria for selection of system gain since the effect 
of root location can be readily seen, and is faster than the rather 
laborious angle summation process necessary in the Evan's Method since in 
most cases only three points must be computed to determine a given curve. 
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IV. Derivation of Transfer Functions of Comronents and Minimum 
Velocity Constant P^auirement 

b .1 Servo A^pl^ fier 

The frequency resDonse of the amplifier was measured and found to be flat 
to + 1 db to 7, COO kc. Since this is beyond the system bandwidth the 
simplification of constant pain will be used throughout this thesis. 

It was found that this pain was anproximately 15C,000. 

Because the input sipnal and the various feedbacks are not summed at a 
common node (see Fip. 11), the pain could not be considered in the 
conventional manner. Rather the pain seen by the various sirnals is 
different. It was found more convenient to consider the pain as a 
function of the variable resistors. This can best be seen oy considerinp 
the equivalent circuit of Fip. 12. 




The pains as a function of the resistors Rj, R^ and R-y is shown in $.1, 

| 

Analysis of the Individual Loop. The entire closed loop eauivalent 
circuit is shown in figure 18 and simplified block diagram equivalent 
is shown in figure 21. 
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h.2 In’efrral Network 

The integral-damping network is a means of compensating the amplifier 
and is shown schematically in Fig. 13. The equivalent circuit is shown 
in Fig. llj. 
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Fig. 1U EQUIVALENT CIRCUIT OF INTEGRAL-DAMP NETWORK 

From Fig. lh it can be shown that the transfer function is: 

Qo _ 5 W.C.i 

~&i ' (R. C t k. ) S * y. ( k. C, + c\ / X, CJ 5 ^ _/ 

k.c. ki C* x, J, Ai c x 

(K-cJ s 



Simol ' fving 



0o _ 

Qi " 



5*-/ ) S + - 

1 ^ <r, /?* <r. / c, X'* <. * 



From Fig. 13 the various values of resistance and capacitance were 

obtained and substituted into the above transfer function which veilds 

/OS 



& — 

when simplified: 

£y *- 



5 *■ /• /z.fz 
/O S 



Factoring: — , , 

(-5 * t* //S^s/ 



A frequency response analysis shown in Fig. 15 revealed that the actual 

transfer function was -^2. _ 2.(,.C$ £ # This transfer function 

< 9.2 <0>*-z){S zo^T 

was used in the analysis. 



It was found later in the analysis that a value of C^ • C^ = ,25^/f 
produced an improved response of the individual loop. Proceeding as 



above the resulting transfer function is 



in this section for convenience. 
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it-3 Servo totor 



The Servo motors Nk 12, Mod 0 used were manufactured by KEARFOTT COMPANY, 
INC., LITTLE FALLS, NEW JERSEY. To obtain the dynamic constants, it was 
recommended that the following relationships be used. 

A r l To#?u£ C Q/V3TAA/7) = ^ CM isajL-r 



k» ( oc/ry c O/V) ^ £* 47+** 
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Using a plot of number of gear teeth versus inertia in 0 2-/*'* the 

total reflected inertia was found and added to the given inertia of the 
motor. The total inertia was found to be J~ ~ X/o * J£-/as*" o< 
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When this transfer function was used in the computation to determine the 
closed loop response and checked against measured response, the results 
were found to be 35 per cent in error. Several different pole-zero con- 
figerations obtained by using various values of tachometer feedback were 
used and the 35 per cent error was found to be consistant. It became 
obvious that the motor time constant was in error. With this in view, an 
investigation to determine the response or a second order system was attempt- 
ed by removing the tachometer and integral networks. Transient resnonse 
curves were obtained at various amplifier gain settings and found to be 
almost identical for any setting of gain. The slight change was accredited 
to the pole of the amplifier at 60C0. From these curves it can be seen 
that ts (settling time) * .355 using the relationship tjjg it was fond that 
o^r /*./ or 2* = 20.2. . This can be seen readily by looking at 

the root locus plot of the second order system shown in figure 16. 
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Fig. 16 ROOT LOCUS OF INDIVIDUAL LOOP AS 
A SECOND ORDER SYSTEM 

From the response 3 

The motor transfer function then becomes: 

/r~ * o 

* J5( 5 +2 8-2) 

As a further check on the accuracy, a frequency response was performed. The 
plot of the magnitude M in decibels is shown in Fig. IT. It was originally 
intended to plot the results on a Nichols chart, transfer these results to 
a Bodie plot to find the break point or time constant. This was found too 
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impractical because "here was enough backlash present to cause the phase 
angle to shift irrationally. Using the value of Hpv from Fir. 17, second 
order curves were utilized to determine / the damping factor. Knowing 
£ and the natural frequency at which it occured, it was found that 

oL - £ n 

o(= .3 7(37.0 - * 

= 27 8 

Averaging the results obtained from the transient and frequency response, 
it was found that the motor transfer function was: 

rr _ ^y- 4 _ 

^ ~ SC Si-23) 

This transfer function was used in all the computations and was found to 
yeild excellent results. 

h .3 Derivation of Theoretical Minimum Kv Requirement 
The velocity constant reauirement was not soecified. Since the system 
under studv is a computation mechanism it was imoeritive that a realistic 
minimum value of accuracv required be determined and used as a oal in 
the analysis. Within reason it would, of course, be desireable to have 
Kv as large as possible in order to reduce steady state errors, but at 
the same time the system need not be designed to have greater accuracy 
than sensing devices within the fire control equipment. 

•f 

In the submarine fire control nroblem, bearing rate in the early stages 
of the approach will generallv ranre' between 0-5 per minute (rarely 



exceeding 2° per minute). As the submarine and target closes, the bearing 
rate will increase, depending upon the relative motion of the two and their 
range. At very close ranges (1000 yards or less) bearing rates may become 
as high as 60 per minute, although this is rare in typical problems en- 
countered by present-day high speed submarines. 

Consider the extreme case which could conceivably be encountered by two 
very-high speed submarines, each travelling $0 knots on opposite courses, 
and passing within 100C yards of one another. Their relative speed is 
100 knots (nautical miles per hour) and their range is 0.5 nautical miles. 
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Individual components within the computer itself have accuracies generally 
of the order of 5 minutes or less. If a steady state error to a velocity 
input of the same mapnitude is specified, then from the relationship: 
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This would represent nearly ideal tracking under the most severe condit- 
ions likely to be incountered in a submarine fire control problem with- 
in the next decade, and is a value toward which it would be desireable 
to work. Obviously, for the great majority of applications K y would 
not have to be nearly as large to obtain 5 minutes or less tracking 
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V. Mathematical Analysis of the System and Comparison with Observed Data 



5.1 Analysis of Individual Loops 

Since the four individual servo loops. Gyro Anple , Enabling Run , Weapon 
Time to Burst , are essentially identical, only one loop will be considered 
in detail. This is justified by measurements made on the other loops. 

For illustration let us consider the Gyro Anple Loop . As previously stated 
in section b.l, the gain cannot be considered in the conventional manner. 
The closed loop equivalent circuit is shown in figure 18 below: 




Fig. 18 CLOSED LOOP EQUIVALENT CIRCUIT OF INDIVIDUAL LOOP 



It is noted that conventional slimming resistors are not necessary. This 
can best be understood by considering the equivalent circuit of the input 
and feedback vernistat in figure 19. 



See Appendix 
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Fig. 19 EQUIVALENT CIRCUIT OF INPUT AND FEEDBACK VERNISTAT 



Further simplifying: 



£ 




Since the reference voltage is variable and R the circuit can 

be arranged in this manner. 
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It is convenient to derive the transfer function of the circuit comprising 
the amplifier and integral damping network, figure 20^ prior to proceeding 
with the analysis. 




Fig. 20 EQUIVALENT CIRCUIT OF AMPLIFIER AND INTEGRAL DAMPING NETWORK 

Breaking the circuit in this manner is justified since neglecting the 
shunting effect of 60K and R^ » 3 megohms results in less error than 
incurred in assumed values of components. 

The resulting block diagram is: 
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Simplifying: 
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Fca = closed loop transfer function of the Amplifier. 
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The Root Locus Eouation is: 
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Tie close loop transfer function of the amplifier network will thpp 

have the form: 4- ~ (±*±S<± tAE^A 

C * {S¥6j 
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The equivalent circuit of the entire loop then becomes: 




The block diagram representation is: 
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By block diagram maniDulation: 
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This simplifies to figure 21* 
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Fig. 21 REDUCED BLOCK DIAGRAM OF INDIVIDUAL LOOP 
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The root locus equation is: 
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Simplifying: 
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The closed loop function vrill have the form: 
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It can be seen from the above equations that varying R2 and Ry varies 
the location of the zero due to the tachometer and simultaneously Ry 
varies the root locus gain* It can be seen that within limits many 
settings of R2 and Ry will give identical results* 



Since Rg and Ry enter into gain as well as zero location, the conventional 
method of considering a eiven pole-zero configeration and varying the 
pain to determine optimum response cannot be used. Further, if a comp- 
osit open loop transfer function is obtained, it can be seen that these 
factors also enter into the commutation of the Velocity Constant, K v . 

N 

Using the relationship F * — the open loop function is: 
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[t is seen that three equations involving the two qualities, response and 
relocity constant, can be easily written. 
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Since the system cannot be optimuraized by varying gain, a plot of these 
three equations will provide a criteria for choosing values of resistance for 
optimum response and maxinaira velocity constant. 

From figures 22 and 23 it can be seen that a near optimra selection of 
resistance is Rx = joo* and Ry - /***" . It is noted that these 
values give Kv in excess of the minimum derived in 1*.3. 
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With these values fixed, consider the loop in three conf igerations : 

1. Single loop - Unity, tachometer, and integral-damp feedback 

2. Single loop - Unity, and tachometer feedback 

3. Single loop - Unity, and integral feedback 



5.2 Single Loc 



Substituting R£ 
becomes: 



Unity, Tachometer, and Intepral-damp FeeriWir 
300 K , R ? - 100 K , and R 6 - 3 Megohms equation ( 5 . 1 ) 
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Equation (5.2), the root locus equation of the amplifier and Integral- 
damp feedback network becomes : 
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The resulting block diagram is: 
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Equation (5.3) becomes: 
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Equatibn (5.U), the final root locus or characteristic equation of the 
loop is: 

y(. 7 (s+i) 
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The root locus plot is shown in figure 2li. The roots are determined from 
K versus the real roots in figure 25. 

The real roots were determined to be S • - 1862 and S ■ - 2.0332 . From 
the characteristic equation, denominator of equation (5.li)> the coeffi- 
cent of the second highest term is 1931 and the constant term is 
6,120,000. Proceeding as outlined in section 3.3s 
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The factored closed loop function then becomes: 
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The response to a unit step input is: 
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The inverse Laplace transform yeilds: 
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Figure 26 shows predicted Q g as a function of time. Figure 27 shows 
the measured response obtained. 



3o2 Single Loop - Unity and Tachometer Feedback 
Removing the integral-damp network, the block diagram is: 




The closed loop transfer function is: 
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Simplified becomes: 
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The root locus equation is: 
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The root locus is shown in figure 28. The Root Locus Gain versus Real 
is shown in figure 29. 
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MEASURED RESPONSE OF INDIVIDUAL LOOPS, TO A STEP INPUT 
USING OPTIMUM VALUES OF INPUT RESISTANCE 
WITH UNITY, INTEGRAL DAMPING AND TACHOMETER FEEDBACK 
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The factored closed loop transfer function is : 
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The response to a unit step is: 
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The inverse Laplace transfonn yields? 
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The predicted output as a function of time is shown in figure 30. The 
measured response is shown in figure 31* 



3.3 Single Loop - Unity and Integral-damp Feedback 
With no tachometer feedback, the block diagram becomes: 




The block diagram simplifies to: 
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The closed loop transfer function then becomes: 
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MEASURED RESPONSE, TO A STEP INPUT, OF THE INDIVIDUAL LOOP 
WITH UNITY AND TACHOMETER FEEDBACK 




Simplified: 
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The root locus equation is: 
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The root locus plot is shown in figure 32. The plot of K versus Real 
Roots is shown in figure 33. The factored closed loop transfer function is: 
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The predicted output as a function of time is shown in figure 3U. The 
measured response is shown in figure 35. 

5.2 Qualitative Stability and Performance Analysis of Overall System 
Section 3*1 outlines the Determinantal Method of analysis used in reducing 
the functional diagram, figure 36 , to the configeration shown in figure 
37. 



It is noted that the Torpedo Lead angle Loop is inherently non-linear 
because of the resolver which is used for a function generator and a 
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MEASURED RESPONSE, TO A UNIT STEP INPUT, OF THE INDIVIDUAL LOOP 
WTTH UNITY AND INTEGRAL- DAMPING FEEDBACK 
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System Functional Diagram 



nulling point. The output from the resolver is Z. X Cos Bb6 - Z / Sin Bb6 
which tends to zero when Bb6 is computed correctly. 

Because of this non-linearity, it was assumed that the inputs Z X Cos Bb6 
and Z T Sin Bb6 were constant inputs and not a function of time. The 
resolver was represented as a summer. Several of the actual constants 
vrere combined with the scale factors of the loop for convenience and simp- 
licity. This is justified since it does not effect the accuracy of the 
stability analysis and it is highly unlikely that quantitative response 
can be predicted with the assumptions made concerning the resolver unless 
a prohibitive number of small perturbations are considered. 

The assumption is also made and justified in section 6.1 that for signals 
which do not saturate the amplifier or motor that the Enabling Run , 

Weapon Time to Burst , and Gyro Lead Angle transfer functions are essentia- 
lly identical. For large signals this assumption is especially invalid 
in the Weapon Time to Burst Loop because of the I4O to 1 gear ratio which 
permits the motor to saturate. 

I 

With these assumptions in mind, the "Standard Block Diagram" becomes 
that shown in figure 37 * Where: 

-A* - & +2. 033*.A5+‘/802.)(S+23+jZ2.(.S‘^5rSS-/ ZI LD 



G. - 

Gj. =? 



k,F, * 




/ 






6Ll . G* 3 - a AG 

O, J / 

G - * 

* 5 / 

Ai . <? 4 - ± Al 



7li 







75 



* 7 - 






/ ^8 




K 

s»1 




4 £ 4l 


*,• 


S /A' 


0£j 


' *>» 


hr = 


Cos 


&£ S 



Lookinp at the outputs 



/fa - s et> 

a Cos S6 
Jri = S/sv &AC 

, Go ( ±J 






* /fc * 



6?*c& - (SV 



/- _ &A& _ Gs~ ^la-f (~JJ 
= ***; ~ 4 



G&c 

~g 7 ^j 



CL+f 



The delta of the system is: 
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can be obtained 






It is convenient to expand the delta so that the A^f 
since it will be used in the numerator of the closed loop function. 
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Simplified 

A - / / Gi G* G 7 t G,C, G<. ^ G, G z Gi G^-Gg ^ 0/ G± GvQ~ 

Ac tf = G, Ci Gj G$ 4 C»Ci Gy. 

The closed loop function becomes 

f a Gf /. C. Cu Oj U g G> t/4. O v J 

/ / a t Cl C, + Cl Cl (7, A Oj-i G, C. Gj Gg t G, G, G y J 

Using the relationship f2 » ^ the open loop function is 

^ 

/r . CV/ G, G i G s Gq y- G , Gi G y j 
/ ^ Gi G y Gj r- Gi Gj Cl 

Since each function G. Gz. etc. has the form of l numerator divided by a 

Substituting and simplifying: 



r> . _ As • 

denominator, let us define — —■ 

8 0*1 



__ / l /, A'i /</y 7 * /(^ ’ a . /^v A A ^7 

<3 />* ^ aT 4 "V /n ^ ^ ^7 ^ 4 44^ 4 4 

In terms of numerators and denominators of the various functions, the closed 
loop function becomes: 



’ _ // ^ /■ /^V y^i /^V ^ 40y 

O.jz o ,/>^ 4- 4 ^7 4? * a / 7 o , 4- 4 ^ ^ >6 4 4* 

/" /V, /t'z, /V> A* 8 Oy 4 ^7 ^ ^ r ^y 
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Substituting the various values for the numerator and denominator of the 
functions, the open and closed loop transfer functions become: 



/ r = (/s, 3oq) v /r Ar x L k\Kti ^ fa 1 ( 5y* ^ oo) ^ 

(S+ IUI)* (5 + 2. 0332J* (S+ S3 * / ZZLS)^ 
[(S^ 33 -J 22. 6 Si' 1 ' J 

r (5 + ItU)*- ( S+ 2.0S3Z) 1 - ( S + 33 +j 22 LD*- 1 
|__V* U3-.300)* 2) *-U3 + Zoo)' X- z ( ^ + I 



p _ (A?; .?«<>;* *T$dU](5/*y* ( j ^A3oo> v 

[/ -S> //U^ 4, ^ S * 2.0 3 31)*- < S + 33 +j Z2.CS)*\ 

■( sy- /Q(,2)* ( 5 V-2. 0332 J 3 " ( Ji 2Z.LS-)*- 
y* uc -joo/*- ts>i/ * ts +2.00)* tr x ( t r + K t rj 
+ (/s; 3oo ) K, k i C Ar, k H - t s+2J+(s+joo)+ 

j m 

Observing equation 5.5 the open loop transfer function is of the form 

F * — [£ *J*+L where n » 0. This indicates, of course, that the 

0 3 h (s +P.K1+P*) 

system is type zero. 



This indicates that the steady state error would be infinite and the sy- 
stem must be compensated in such a way to make it a type one if ramp 
inputs of large magnitude or long duration are expected to be incountered. 



/<v = 5 s. 



J -f o 




IS+£,J 

(S + P.U&+/U • • 




The steady state error to a ramp is: 




CO 



& = S3. OPS OS AtAAlP 
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Since, of necessity, certain linerizations were made; therefore, no 

conclusions can sensibly be made concerninr stability except that the 

necessary condition that all the coefficients of the characteristic 

equation are positive is met. This being the circumstance, a frequency 

response analysis was made. The results were plotted on a Nichol's 

chart. It was found that there was a gain margin of 26 db and a phase 
o 

margin of 72 . This, of course, shows that the overall system is highly 
stable. The results obtained from the Nichol's chart were then plotted 
on a Bodie diapram, figure 38 » to determine the dominant time constant 

of the system and to acertain experimentally that the system was actually 

^ . 

type zero. These results are discussed in section 5.3-b. 

5.3 Discussion of Analysis and Observed Data 
5.3a Individual Loop 

The individual loop was analyzed in three possible configerations, unity, 
tachometer, and integral-damp feedback; unity and tachometer feedback; 
and unity and integral-damp feedback. 

By using figure 22 and 23 it was found that near optimum values of 
resistance were Rj ■ 300 K and Ry * 100 K. It was found that the value 
of R^ used had little effect on the response of the system; however, 
considerations of noise in the system dictated a value of approximately 
three megohms. Using these values of resistance, computations were made 
to predict the response of the loop. The predicted and measured response 
is shown in figure 26 and 27 respectively. 

Prior to the determination of the values of resistances, the best response 
obtained was a critically damped case with a settling time in excess of 
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.21 i seconds. It can be seen that the predicted response and measured 
response agree exactly and a settling time of .1$ seconds was obtained. 

This shows an improvement in response of 37*5 per cent. 

The analysis of the loop without the integral-damp network was performed 
in section 3,2. The measured response is shown in figure 31. Comparing 
this with the predicted response in figure 30, it is again seen that the 
measured and predicted values agree. It is also noted that an improvement 
in response was obtained with a settling time of .11 seconds. This, of 
course, is desireable; however, at steady state there was audible noise 
caused by gear chatter. This noise did not show up in the output but 
would cause excessive wear in the gear train near the motor j therefore, 
the integral-damp network is considered desirable. It was found analy- 
tically that by changing capacitors Cj and (See figures 13 or lii) to 
.2$ jurf a predicted and measured response with a settling time of .12 
seconds was obtained without any associated noise problems. The theo- 
retical pole location is calculated in section 1.2. 

The last configeration considered involved unity and integral-damp feed- 
back. From figures 3b and 35 it can be seen that the predicted and 
measured response are again in agreement except the peak overshoot is 
somewhat greater in the actual system. Settling time is .3 seconds. 

It is noted that substantial tachometer feedback must be used since for 
smooth computation each loop must be critically damped. 

5.3b Overall System 

As discussed in section 5*2, the overall system is highly stable with a 

o 

phase margin of 72 and a gain margin of 2Ij db. The data from the Nichol's 
plot is used to construct figure 38* From this it is noted that the 
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first asvmptote has a zero db slope; therefore, type zero system. The 
time constant agrees reasonably with the actual obtained from the meas- 
ured response shown in figure 39. “Hie theoretical settling time is .6 
seconds and the measured is .7 seconds (See figure 39). The response 
to a ramp input is shown in figure bO. 

% 

Since the theoretical computations and frequency response indicated that 
the system was type zero with a finite error to a step input and an 
infinite error to a ramp input, it was considered advantageous to meas- 
ure these errors. The inputs and resulting errors are shown in figures 
bl ancl U 2. 

The results are therefore conclusive that the system is in fact type zero. 
While the multi-loop system is unmodified (i.e., a type 0 servomechanism), 
there will be a constant error to a step and a corresponding, increasing 
error to a ramp input. The measurements in figures hi and b2 show the 
magnitude of the errors involved. The first is a measurement of the error 
versus time for the step and the second is the measurement of the error 
versus time corresponding to the ramp. From the first the error ratio 
to the step input can be calculated to be: 

- . oos~ 

7. 5" rai. 

( 

In a typical submarine problem, the total bearing change from initial 
contact to the firing point is approximately 1$Q degrees. This means 
that, under these circumstances, the error at firing point will be: 

(-OOS •) - . Q° 

One degree ■ 17.6 yards/1000 yards. Therefore, the following errors 
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MEASURED RESPONSE OF THE OVERALL SYSTEM TO A STEP INPUT 
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MEASURED RESPONSE OF THE OVERALL SYSTEM TO A RAMP INPUT 




MEASURED ERROR OF THE OVERALL SYSTEM TO A STEP INPUT 




MEASURED ERROR OF THE OVERALL SYSTEM TO A RAMP INPUT 



woi Id result with a computed error in the Gyro Lead Angle Computation 
Loop of 0.8°: 

at: 1,000 yds. : lii.05 yds. 

2.000 yds. : 28.10 yds. 

3.000 yds. : 1*2.15 yds. 
i*,000 yds. : 56.20 yds. 

These errors are, of course, in addition to any other errors generated 
in other parts of the computer. It can be readily seen that if fire 
control problems of short range with the resultant large ramp or long 
ranges where tracking is done for long periods of time, a significant error 
may develope. 

Looking at equation (5.5)>the open loop transfer function of the system, 
it is seen that an integration is required or an n S" is needed in the 
denominator to have a type one system. It appears that a feed forward 



loop added to the block diagram, figure 37, would accomplish the desired 
results. The resultant standard block diagram is shown in figure 1*3. 




Fig. 1*3 STANDARD BLOCK DIAGRAM WITH COMPENSATION 
NECESSARY TO PRODUCE A TYPE ONE SYSTEM 
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The determinant for delta becomes 
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The determinant reduces to 

A = // dadvd, 74 da G, dl 74 C, / <?, d, < 7 , dy 74 < 7 , C, d„ dy 

y <?' <da <?v dr y <£ C C C C* 

Aa^ - Q <j t da dv Oy dy 74 da da dV dy 14 d da dv dp r d da da dj-Cy 



The closed loop function is 

-- jGt / da da <?£ dy / da dyd, dV ^d,da dv <?,■ / dl da da <?, d / ) 
/d 74 da dv dV 74 da da dy 74 ^ 74 da dv dy 74 da dy da dy 7 

[r 4 da da dv da- 74 d, da d, dp d, J 



In terms of numerator and denominators 



/I'V A Pi. Pi P« 4^7 ^ ^ A/V -&7 My P/ Pj Py Pc P)q 

Mx. y^i Me M*; P* Py Py P? <0$ t 4 - M/ M*. My Asj- Pi Pr PyP$ 0 9 

+ M, M A/j MrA® P± Pc Ph_ 

Pi P* Pv Py Pq P 7 Ps P? ^ M*. M A, p> Jj- A- P 7 Pe 
A Mi. A'j Me />, P« Pr P7 P$ Py * ^ P* Pi Pj Py Pj- Pr P? P9 
+ Mi AV *7 Pi Pj Py Pe P$ * Ai AV M r Pi Py Pj-Pt 
/ /V/ Mi /t'v Asy Pj Pc Py Pq Pf A AlS* AfM.A'sMs'AijfP^P^Pj/^ 
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The open loop function is 



A? 0, Ox O 3 Ay Os' At A*7 tS$ ^ A*. Ay A> Ay A/ Ai Ay- At Ay 

/ A* A> A<. Ay A/ Ay Aj~ A 7 Aj 7*- A/ A*. Ay Ar Ay A«, Ay Ay Ar 

7*" A/ A'*. A'j Ay* Ay Ay At A 7 Ay 

~ 0 > Ox Oj Ay Aj - At A* A/ A* * A'*. Ay A> A, Ay Ay- At A, A, ^ 
/■ Ai Ay At A^ Ay Ay- A 7 Ajp A7 

It can be seen that the necessary and sufficient condition for a type one 
system has been met since 0/y * S >*iZ J can be factored from each term. 

Where G ^ • *2 is of the form: - s( s+</ • 
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VI. Conclusions and Recommendations 



6.1 Single Loop 

The single loop designed and constructed by Librascope Division of General 
Precision Corporation will produce satisfactory response without further 
changes in design. However, by judicious selection of the various input 
resistors, an improvement is possible of 37.5 percent using settling 
time as a criteria. It can be seen from the root locus, figure 2h and 
K versus Real roots, figure 25 that the system is absolutely stable for 
all values of gain. It can also be seen that the gain can vary appreci- 
ably without effecting the response. Should the gain become very low, 
an appreciable residue will be associated with the root near the zero 
at S * -2. Tachometer feedback is essential to provide a critically 
damped system for computation. Integral-damp feedback slows down the 
response by 26.5 per cent from that obtained without it; however, the 
damping is required to prevent high frequency oscillations in the gear 
train during steady state operation. In selecting values of capacitance 
(C 3 ■ ■ .25 M’f ) in figure lb, the deterraental effect on response 

due to the integral-damp feedback is eliminated. The recommended changes 
improves the response a total of 5 b . 2 per cent over the best obtained 
by trial and error on the physical system. 

6.2 Overall System 

The overall system is type zero as designed. This is satisfactory for 
most applications intended since the error involved is small; however, 
should large ramp inputs or tracking problems of long duration be encount- 
ered, it is possible to generate excessively large errors. The solution 
to this problem is a modification involving a motor and amplifier in a 
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feed forward loop. This loop will be physically connected at the error 
input to the Torpedo Lead Angle and feed forward to the output of Bb6. 

This is represented in the modified standard block diagram, figure L?. 

An alternate method which may be prefered, since the primary function 
of the loop is to mechanize an equation, is to determine a mean ramp, 
calculate the various outputs as a function of time, and scale the various 
loops to reduce or eliminate the resultant error. A check must, of course, 
be made over the full range of ramps expected, after the scaling has been 
accomplished, to insure that the error does not become excessive else- 
where. ‘ 

Measured frequency and transient response indicates that the system operation 
will be satisfactory with .7 seconds settling time to a unit step, gain 
margin of 2h db, and phase margin of 72°. 
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VII. Recommendations for Further Investigation 






Should it be deemed desirable to predict accurately the response to a 
complex multi-loop system utilizing resolvers, it would be advantageous 
to use the linearizations made herein, consider small pertubations, and 
determine their accuracies obtained bv measurements on the physical system. 
This would provide a criteria for evaluating future analysis which may be 
conducted during design work on other systems. 

In addition it may be advantageous to consider linearization technicues 
other than the ones used herein. 

If it is reouired that the overall svstem be optimized, an analogue computer 
study is warranted. 
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APPENDIX I 



Induction Resolver 

An induction resolver is an electromagnetic device employed as 
a basic element in analog computers, automatic control systems, 
angle data transmission systems and plan position indicating 
radars. It is used for performing trigonometric computations. 

The induction resolver is essentially a transformer with rotaiy 
variable coupling between primary and secondary windings (Figure 1). 
The windings are distributed in the various stator and rotor slots 
in such manner as to obtain a true sinusoidal relationship between 
output voltage and angle of rotation of the rotor. Since the 
units are used at audio and power frequencies, it is customary 
to use a high-permeability laminated core similar to that used in 
an induction motor. 




It is customary to treat the stator winding as the primary or 
input winding, and the rotor winding as the secondary or output 
winding. If, as in Figure 2, there are two secondary windings 
at right angles to each other (in "space quadrature"), one sec- 
ondary will have induced in it a voltage which is proportional 
to the sine of the angle of rotation between it and the primary, 
and the other secondary will have induced in it a voltage which 
is proportional to the cosine of the angle of rotation between 
it and the primary. 



To make induction resolvers even more useful, two windings are 
provided on both the stator and rotor. The pairs of stator wind- 
ings and rotor windings are located at right angles to each 
other as shown in Figure 3. This permits computation of cosines 
as well as sines, and even more complex problems involving both 
sines and cosines can be solved when both stators are excited. 

In addition to the resolution of input voltages into sine and 
cosine components, resolvers can perform the synthesis of resul- 
tants from components and the rotation of coordinate systems. 



\ 

\ 




//V/’l/ T 3 QU7~ 7~ 3 



V 0 

Q C Q 



Fig. 2 - Resolver with two Rotor Windings 
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Fig. 3 - Resolver Input and Output Relationship 

As shown in Figure b, currents in the two stator windings 
set up magnetic fields at right'-angles to each otter to form a 
resultant flux vector. This flux vector may assume ary angular 
position, depending upon the amplitudes and polarities of the 
input currents. Its magnitude and proportional to the vector 
sum of the instantaneous values of the input currents, and its 
angular position is equal to the arctangent of the ratio of 
the two currents. For AC input voltages within the operating 
frequency range of the resolver, the rotor voltages are propor- 
tional to the magnitude of the flux vector and the sine or cosine 
of the angle of rotation of the rotor winding with respect to 
the vector. The output voltages may be expressed in terms of 
the input voltages and the rotor angular position, by the fol- 
lowing equations: 
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( 1 ) 

( 2 ) 



F# - F Coj 9 - £~ 3 3//v & 

F# = F CojQ + FI 5/sro 

‘ *■ J/ J 2 

The subscripts refer to the various rotor and stator windings 



and O is the rotor angular position. 

— Xt. 
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Fig. h - Stator Flux Vectors, 



While other types of operations are possible with a resolver, 
its use has been principally with three types of inputs, namely; 
(l) sine wave voXtages of the same electrical phase, (2) sine 
wave voltages shifted 90 degrees apart, and (3) irregular wave 
shapes, principally of the sawtooth or square wave variety. For 
the first type of operation, the output voltages remain fixed 
in electrical phase shift and their magnitudes vary in accordance 
with equations (1) and (2). For the second and third types of 
operation, the output voltages may also be determined by equations 
(1) and (2) if the instantaneous values of the voltages are taken 
into consideration. For the case of the resolver with two inputs 
90 degrees apart in electrical phase, each of the two rotor output 
voltages remains constant in amplitude, independent of rotor an- 
gular position. However the phase of the voltages varies contin- 
uously with rotor angular position. The phase shift between the 
two output voltages remains fixed at 90 degrees. In this mode of 
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operation, then, the equations for the output voltages are: 







L*. 


(?; 


r*. - 




/o +?o° 


(M 



For the third case, where irregular wave shapes are used, op- 
eration most frequently falls into two categories. In phase, 
square waves or flat topped sine waves are used in some computing 
applications. The performance is adequately described by equations 
(1) and (2) without modification if the two input wave shapes 
are the sane, which is usually the case. 

Sawtooth wave shapes are frequently used in radar sweep appli- 
cations. In the simplest case, a sweep voltage is applied to 
one stator and the two output voltages from the rotor windings 
are the same sweep voltage but with an amplitude proportional 
to the sine or cosine of the rotor angular position. If these 
two voltages are applied to the "X" and "Y" deflection plates, 
or coils of a cathode ray tube, they cause a radial trace to be 
formed. The trace will assume the same angular position as the 
resolver rotor. This type of operation finds application in 
plan position indicator types of radar displays. 

While the arrangement shown in Figure 5a is suitable for schema- 
tic representation of a resolver, tie arrangement shown in Figure 
5b is more suited for block diagrams and functional schematics. 
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The diagram, equivalent schematic and equations are shown for 
clockwise rotation (facing the shaft end) and apply to all 
standard types of resolvers. Vhere it is desired to change the 
direction of rotation of the resolver, it is only necessary to 
interchange input A with B and output C with D. This can be 
noted n n the block diagram by a negative sign preceding e . 




Figures 6a and b. One stator winding is excited with a voltage 
proportional to the hypotenuse of a right triangle and the rotor 
is positioned to the angle <3* . The two output voltages are pro- 
portional to the two legs of the triangle. 



T/Q*<j£' 7 ~ 




Fig. 6 - Elementary Resolver Application 
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In Figure 7 the problem is similar, but the hypotenuse and one 
of the legs, Y are known. A servo amplifier and motor are used 



t 9 position the resolver rotor until the voltage, is del- 

ivered by the resolver equals the known voltage Y. The voltage 
is then proportional to the length of the other leg of the 
triangle, and the angle © is the angle between the second leg 
and the hypotenuse. A solution of this type may be used in de- 
termining elevation angle and ground range when slant range and 



depth are known. 



7 *3 » * * 




^ ‘ ^2 I © ~ 




Fig. 7 - Angle Confutation 



In Figure 8 , the inputs are two vectors 90 degrees apart. TJie 
servo positions the resolver rotor until one output is zero. 

At this position the other rotor being 90 degrees away, has maxi- 
mum coupling with the flux vector formed by the X and Y inputs. 
Its output is therefore proportional to the vector sum of these 
two inputs. The rotor angular position is the angle between one 
leg and the hypotenuse of the triangle. 





Fig. 8 - Vector Addition 
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In the rotation of coordinates, both stators and rotors are 
used. With rectangular coordinates X and Y as inputs to a resol- 
ver, as in Figure 9, the outputs are also rectangular coordinates 
of the same point referenced to a new set of axes rotated with 
respect to the first. The operation if fully described by the 
equations in Figure 3 which will be recognized as those required 
for coordinate rotation. 

The actual operation can best be understood by breaking it into 
two steps. The X-Y coordinates impressed on the two stators 
form a flux vector in the resolver. This vector is independent 
of rotor position. Its angular position is the arctangent of Y/X. 
Each of the rotor windings develops a voltage proportional to 
the sine or cosine of its angular position with respect to the 
flux vector (not with respect to zero). If the reference is made 
with respect to zero then the outputs follow the equations 
in Figure 9. "Tie rotor voltages are the X 1 , Y* components of 
the flux vector. The angle between the X', Y* coordinate ref- > 
erence of the inputs is the sane as the rotor angular position. 



V. 
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Fig. 9 - Coordinate Rotation 
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INDUCTION POTENTIOMETER 



An induction potentiometer, or linear synchro transmitter, may 
be considered a special kind of resolver which provides 
accurate linear indication of shaft rotation about a reference 
position in the form of a polarized voltage whose magnitude is 
proportional to angular displacement, and whose phase relation- 
ship indicates direction of shaft rotation. The principal differ- 
ence between a resolver and an induction potentiometer is that 
the latter's output voltage varies directly as an angle and not as 
the sine of that angle. Its rotation is usually limited to less 
than +90°. 




Fig. 1 - Induction Potentiometer Voltage Output 
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These devices are analogous to resistance potentiometers, but 
since they are induction type components they have less re- 
straining forces acting upon their rotors, and hence are capable 
of providing better resolution. Naturally, induction poten- 
tiometers do not require sliders such as those used in resistance 
type 8. Therefore, circuit interruptions are eliminated, no wear 
occurs as a result of rubbing parts, and accuracy is consequently 
continously maintained at the original level throughout the 
operational life of these types of components. These advan- 
tageous features are virtually a necessity in certain applications. 
For example, in gyroscope systems where low restraining torques 
and low pick -off angular errors are required, wide use is made 
of these devices. In general, induction potentiometers find use 
in applications where resistive potentiometers are impractical, 
principally because of the following features: (l) induction 

potentiometers, having no wiping contacts, may be used as 
gyroscope pick-offs since they contribute less spurious friction 
torque: (2) input and output are isolated: (3) resolution is 

infinite: (U) noise level is low: and (5) the total angle of travel 

is limited to less than 180° 

Constructed very much like a resolver, induction potentio- 
meters differ in that their windings and slots are not uniformly 
distributed, but instead are deliberately modified to produce a 
linear output. Normally they have only one excitation winding 
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and one output winding, the former usually carried on the rotor 
and connected by means of two sliprings and two brushes. 



The use of induction potentiometers as linear pick -offs may 
be considered as a special resolver application in which mutual 
coupling varies not as *0, but as 2 Q . With zero input im- 

pedance and infinite output impedance, the voltage gradient. 



H, may be defined as H * — . 

fjQ 

fer function may be expressed as: 

jf: _ 

^ / 

/V 



Under this condition the trans- 



- HO 



is -H x 2r 



where 4 * secondary impedance 

* primary impedance 

* load impedance 
i± « input impedance 
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